Fluctuations in classical sum rules 



O 



Q 

U 
d 



> 

in 
cn 

in 
o 

On 
O 



X 



John R. Elton 1,2 , Aral Lakshminarayan 1,3 , and Steven Tomsovic 1,3 ' 4 
1 Max- Planck- Institut fur Physik komplexer Systeme, D-01187 Dresden, Germany 
2 Department of Physics, Cornell University, Ithaca, NY 14583 USA (current address) 
3 Department of Physics, Indian Institute of Technology Madras, 
Chennai, 600 036, India (permanent address) and 
4 Department of Physics and Astronomy; Washington State University; Pullman, WA 99164-2814 (permanent address) 

(Dated: July 2, 2010) 

Classical sum rules arise in a wide variety of physical contexts. Asymptotic expressions have been 
derived for many of these sum rules in the limit of long orbital period (or large action). Although 
sum rule convergence may well be exponentially rapid for chaotic systems in a global sense with 
time, individual contributions to the sums may fluctuate with a width which diverges in time. Our 
interest is in the global convergence of sum rules as well as their local fluctuations. It turns out 
that a simple version of a lazy baker map gives an ideal system in which classical sum rules, their 
corrections, and their fluctuations can be worked out analytically. This is worked out in detail for 
the Hannay-Ozorio sum rule. In this particular case the rate of convergence of the sum rule is 
found to be governed by the Pollicott-Ruelle resonances, and both local and global boundaries for 
which the sum rule may converge are given. In addition, the width of the fluctuations is considered 
and worked out analytically, and it is shown to have an interesting dependence on the location of 
the region over which the sum rule is applied. It is also found that as the region of application is 
decreased in size the fluctuations grow. This suggests a way of controlling the length scale of the 
fluctuations by considering a time dependent phase space volume, which for the lazy baker map 
decreases exponentially rapidly with time. 

PACS numbers: 05.45.-a, 05.45.Ac,45.05.+x,03.65.Sq 



I. INTRODUCTION 

Classical sum rales play an important role in a number 
of physical contexts. It is interesting to have of a 

non-uniformly spreading chaotic system where their lim- 
its, fluctuations, and convergences can be studied analyt- 
ically, which is given in this paper. We focus on the sum 
rule of Hannay and Ozorio dc Almeida, who were moti- 
vated by a desire to understand the two-point quantum 
density of states correlator [l[. Their derivation relied 
on the Principle of Uniformity, which states that the pe- 
riodic orbits, weighted naturally arc uniformly dense in 
phase space. There are many others; see, for example, the 
sum rules in Refs. 043 which involve return probability, 
the connection of two points in coordinate space, and 
fixed orientations of initial and final velocities, respec- 
tively. They arise in mesoscopic conductance, diffraction 
contributions to spectral fluctuations, and chaotic quan- 
tum transport, also respectively. 

In one possible form of the Hannay-Ozorio sum rule @ 
for chaotic systems, it is the inverses of the stability 
matrix determinants for the periodic orbits of period t, 
|Det (Mr - 1)| , which are summed. For completely 
chaotic systems the finite-time stability exponents, which 
roughly determine the value of each determinant, con- 
verge as t — > oo to a single set. Nevertheless, the de- 
terminant's fluctuations from one periodic orbit to the 
next grow without bound in the same limit. This feature 
is made even more curious by the expectation that the 
Hannay-Ozorio sum rule converges exponentially rapidly. 
See a preprint by Pollicott @ for a theorem regarding the 
sum rule's convergence. 



It is known that smooth classical functions have ex- 
ponential decays in fully chaotic systems toward their 
ergodic averages, which are governed by the Pollicott- 
Ruelle (PR) resonances |7Hl0l] associated with the 
Perron-Frobenius operator. It is natural to ask whether 
the kinds of classical sum rules encountered always con- 
verge to their limiting values with corrections that decay 
exponentially according to the leading PR resonances. 
The work of Andreev et al. [ll| would suggest that as long 
as there is a gap in the spectra of the Perron-Frobenius 
operators, there is no other possibility. 

A first step is taken here in approaching this general 
line of questioning by considering a form of the Hannay- 
Ozorio sum rule for fixed period in maps. Its convergence 
rate and local fluctuations are considered in detail in a 
simple version of a lazy baker map introduced by Balazs 
and one of us ( AL) [l2|, [HI • Although the usual baker 
map is especially simple, it is also non-generic in that 
its lacks the essential fluctuations of interest, whereas 
the lazy baker map possesses fluctuations and still turns 
out to be analytically tractable. The fluctuations of 
the inverse determinant were studied briefly in Ref. flij ]. 
where they were demonstrated, not surprisingly, to be 
extremely sensitive to even tiny islands of stability. In 
the form of the lazy baker map studied here, all orbits, 
with one exception of a marginal orbit of period two, are 
unstable, and there is no ambiguity concerning whether 
the system is completely chaotic. 

The organization of this paper is as follows: in Sect. |TT] 
notations are specified and definitions of various quanti- 
ties of interest are given. Section IIIII introduces a simple 
version of a lazy baker map (named the SR map), which 
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is studied in detail, and gives basic counting results for 
fixed points and a method of subdividing the phase space 
into local regions. This is followed in Sect. lIVI bv a deriva- 
tion of the main results for the local fluctuation and con- 
vergence properties of the Hannay-Ozorio sum rule. 



II. MEASUREMENTS OF CHAOTIC SYSTEMS 

Dynamical systems theory has generated a number of 
ways to specify the complexity of a chaotic mapping. 
Three of the more familiar concepts to physicists are the 
topological entropy, hx, the metric or Kolmogorov- Sinai 
(KS) entropy, and the Lyapunov exponent, Xl flH IT3|. 
The topological entropy is designed to measure the infor- 
mation content of the optimal partition of the dynamics. 
It turns out for that a class of systems known as Axiom 
A, the limiting value of the exponential rate of increase in 
the number of fixed points with iteration number gives 
this entropy. The KS entropy can be thought of in a 
similar way, except that it is weighted. And finally, the 
Lyapunov exponent measures the exponential separation 
of neighboring initial conditions. 

For the purpose of studying the fluctuations of clas- 
sical sum rules, the main quantities of interest tend to 
be the number of fixed points, their finite-time stability 
exponents, and their probability densities and moments, 
all of which can be considered in both a global and local 
phase space context. We do not worry as to what exact 
relations exist between each of these measures for a given 
system, but they are closely related where they do not 
have an identical counterpart and it is useful to relate our 
results to some of these quantities when they are known. 



A. Basic quantities 

The notation N T denotes the number of fixed points 
at integer time r taken over the full phase space. This is 
distinguished from a local count of fixed points by writing 
N T (s, k) where the parameters s and k conveniently spec- 
ify the location and size respectively of the local phase 
space volume in question for the SR map. Additionally, 
of great importance in this work are the sums over fixed 
points appearing in one form of the Hannay-Ozorio sum 
rule, which are given the notation 
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(global) 
(local) (1) 



where "f.p." denotes fixed points in the specified region 
of phase space and 



is the Jacobian stability matrix along the trajectory fixed 
by the set of iterates {qi,pi} of the initial conditions 
(qo,Po); the notation for initial conditions is mostly left 
suppressed. A related quantity which is sometimes of in- 
terest is the finite-time stability exponent, given here for 
a map with single position and single momentum coordi- 
nates 
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where tr(...) denotes the trace operation. For long 
times, the approximate relations for A tend exponentially 
quickly to the first relation. 



B. Convergence and measuring sum rule 
fluctuations 



The Hannay-Ozorio sum rule in this context and no- 
tation (also referred to as the uniformity principle) reads 
F T (s, k) — > Vk where Vk is the phase space volume over 
which the fixed points are summed. The absence of an s 
dependence is the absence of a dependence on the loca- 
tion in phase space. This result, which holds in both the 
local and global cases, emerges in the limit of long times. 
The simplest determination of its convergence amounts 
to calculating the leading corrections to the sum as a 
function of {r, s, k}. The expectation is that it should 
decrease exponentially with t if {s, k} are held fixed, al- 
though it is not obvious to us, a priori, at precisely which 
rate and with what length oscillations. Ahead, it turns 
out to be governed by the leading Pollicott-Ruelle reso- 
nance, which interestingly enough, has a real part equal 
to the Lyapunov exponent in the SR map. 



1. Local convergence boundary 

It is also possible to consider the convergence with in- 
creasingly smaller local regions in the phase space. By 
controlling the size of a region, the number of terms con- 
tributing to the local sum can be tuned for a given r. 
Given that the individual stability determinants fluctu- 
ate with ever increasing width as r increases, the question 
becomes, "at which time on average do the corrections 
become smaller than the sum rule expectation, i.e. the 
local phase space volume?" As the regions shrink in size 
this time extends later, thus making it possible to find a 
shrinking volume as a function of r that offsets the ex- 
ponentially decaying corrections. One can think of this 
"V(r)" as the boundary for the local application of the 
sum rule to be converged. Both the global and the local 
boundary of convergence is given ahead. 



2. Moments 

For the subregions of phase space, subtracting the local 
volume from the sum itself, call it F T (s, k), gives the lead- 
ing corrections and fluctuating components of the sum 
rule. A probability density for the values it takes on at 
fixed time over all regions can be defined, Pp [x), which 
carries all information about convergence, local or global, 
and fluctuations. In particular, our focus in this paper 
is on central moments of the density. A very important 
case is the mean square deviation 

a 2 (F T ,k) = (F T (s,k) 2 ) s = Jdxx 2 Pp T (x) (4) 

where x — x — x and for which an asymptotic expression 
is derived in the case of the SR map. More generally, the 
n th central moment is 

M n (F T , k) = (F T (s, k) n ) s = J dx x n P K (x). (5) 

The moments 

M n {e x \k) = (e"M'.*^) s = J dx x"P cxp(AT) (x) (6) 

which by Eq. §S§ is associated with the probability den- 
sity of the finite-time stability determinants, are dis- 
tinctly different from the moments of the sum rule fluc- 
tuations. Both cases are treated in this paper. However, 
as also shown, the probability density P F (x) asymptot- 
ically tends to a Gaussian density and only the first two 
moments (cumulants) are considered in detail. 
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FIG. 1: Partition of the unit square into four subregions, TZi 
through IZ4, (upper panel) and where each subregion maps 
after one iteration (lower panel). As shown, IZi and IZ2 are 
squeezed and stretched whereas 7^.3 and IZs, are rotated. 



III. THE LAZY BAKER SR MAP 

The Hannay-Ozorio sum rule fluctuations may be 
worked out exactly for the case of a simple dynamical 
system which is a modification of the usual baker's map. 
Lazy baker maps have previously been introduced [T2[ 
as a class of 2D area- preserving maps. We study here a 
particular case called the SR map (stretch-rotate) which 
is chaotic over the whole measure and is defined on the 
unit square in the usual position, momentum coordinates 
as follows: 



q' 


= 2q \ 


p' 


= P/2 J 


</ 


= 1-P 


p' 


= q 



if 0<q<I/2, 
if 1/2 <q<l. 



(7) 



The action of the map can be pictured most easily by 
splitting the unit square into four equal subsquares, IZi 
through 7^4, shown in Fig. [T] The region TZ^ is rotated 
uniquely into IZ3 on the next iteration and region IZ3 is 
rotated into 7Z2- On the left half of the square, points in 
IZ2 and IZi are compressed by a factor two along the p 
axis and stretched by the same factor along the q axis. 



The combination of rotation and stretching in the same 
dynamical system gives rise to the possibility of nonuni- 
form hyperbolic motion or even non-hyperbolic motion. 
For the SR map defined as above with the vertical cut in 
the middle of the square at q = 0.5, the motion is of the 
former kind. As the cut is moved to the right it ceases 
to be completely hyperbolic at the golden mean. In this 
paper we will only study the case when there are equal 
regions stretching and rotating, which is arguably the 
simplest "exactly solvable" model of nonuniform hypcr- 
bolicity in an area preserving map. It admits a Markov 
partition of phase space and the dynamics is one of sub- 
shift of finite type on three symbols as shown in [l3[ . The 
atoms of the partition are A = 1Z± B = IZ3, and 

C = K 4 . 

The transition matrix is 




(8) 



whose ij element, tij, is I if there is a transition from i 
to j and otherwise. Here i,j € {A, B, C}. This has 
only topological information. Let be the fraction of 
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atom i in atom j, on one evolution of the map. This gives 
the transition probabilities for the 3-state Markov chain 
that the SR map is equivalent to. The transition prob- 
abilities are p A A = 1/2, vac = 1/2, Pba = 1, Pcb = 1, 
the rest being zero. The Markov matrix is then T\ whose 
matrix elements are tijPij while for fluctuations to be 
studied below the matrix Ti whose elements are tijpfj is 
also useful. These matrices can be used to study the uni- 
formity principle at the global as well as the local scales 
as is shown below. To study what happens on restriction 
to smaller areas whose size can be controlled, as well as 
to find the actual locations of the orbits, it is useful to 
use a binary representation given ahead. This is not a 
symbolic dynamics in the sense that the dynamics is no 
more a left shift. However the dynamics is an eventual 
left shift even in the binary representation and is used 
extensively below. 

Any point in phase space can be represented as a bi- 
infinite binary string p.q representing its position q and 
momentum p. The binary string representing the first m 
bits of the position coordinate is labeled 7 m . The quan- 
tity to also corresponds to the number of times an orbit 
visits the stretching region (left half) of the square after 
some number of iterations. The rules for the mapping 
equations on the binary string are given in fl2| and sum- 
marized here: if the most significant bit of position is a 
0, the dynamics is that of a left shift. If the most sig- 
nificant bit is a 1, position and momentum coordinates 
are interchanged and in the new momentum coordinate 
all 0's and l's are switched. For example, consider the 
period 3 orbit starting at (qo,po) = (2/3, 1/3): 

(2/3, 1/3) -> (2/3, 2/3) -> (1/3, 2/3) -> (2/3, 1/3). 

The binary representations for the starting coordinates 
are q = .10, p = .01, where the underline indicates infinite 
repetition, and under the dynamics this point maps as 

OLIO -> 01.01 -> 10.10 -» 01.10 (9) 

while under the symbolic dynamics this orbit is CBA — > 
BAC -> ACB . 

Making use of the symbolic dynamics, it is possible to 
prove that the map contains a dense set of periodic or- 
bits and is hence an ergodic transformation. It is also 
a straightforward argument to see that the positive Lya- 
punov exponent Al for the SR map is 

A L = ^ln2. (10) 

This may be seen as follows: The unique smooth invari- 
ant density is uniform on the phase space. Therefore 
crgodicity implies that at a large number of iterations 
a typical orbit spends equal amounts of time on the left 
and right halves of the unit square. Since points along an 
orbit which are on the left half arc stretched by a factor of 
two and points on the right half are not stretched at all, 
the Lyapunov exponent will be the average of In 2 and 0. 
The known symbolic dynamics, or binary representation, 
also allow the enumeration of all periodic orbits of any 
period as seen in Sect. MI Bl 



A. Periodic Orbits and Stability 

To begin the discussion of the periodic orbits, first note 
that there exist two exceptional periodic orbits on the 
boundary of the square: a period 1 fixed point at the ori- 
gin and a period 2 orbit between the points (1, 1/2) and 
(1/2, 1). All other orbits must pass through the interior 
of the rotating region IZn and it is thus sufficient to count 
orbits originating in TZ&. The orbits come in two types 
depending on the parameter henceforth called j, which 
is the number of changes from to 1 or 1 to in the 
binary string 7 m representing the first m bits of the posi- 
tion coordinate. If j is odd, the orbit may be represented 
as p.q = r y m -"fm- If j is even the orbit may be written as 
P-1 = 7m7 m -7m7 m where 7 y m denotes the complement of 
7 m . From the rules given for the symbolic dynamics one 
finds the period r of an orbit in terms of its 7 m string to 
be 

T = 2j+TO + 2 (11) 

where j may take any integer value from to m — 1. 
Reference [12[ may be consulted for more details. 

It is also useful to realize that it is possible to 
translate from the binary to the symbolic dynamics 
and vice-versa. Sticking to orbits originating from 7?_4 
these are of the form ••• 0. IX2X3 ■■■ . This translates 
by replacing every transition (which is a 0-1 or 1-0 
"bond" in the binary representation) including the ini- 
tial 0.1 by CBA and every other type (i.e. 0-0 and 
1-1 bonds) by AA. Thus for example the orbit with 
the binary representation ••• 0.10100 •• • translates to 
• ■ ■ CB ACB ACB ACB AAA 

The Jacobian stability matrix for a single time step 
is dependent on whether the point in question is in the 
left or right half of the unit square. Denoting Mj, as 
the stability matrix for points in the left half and M# as 
the stability matrix for points in the right half we have, 
following from the definition of the map, that 

M L = ( 2 Q ^ , M R = - 1 ) . (12) 

The stability matrix along a trajectory Eq. ^ may then 
be calculated explicitly as a product of the matrices 
Eq. (|12p . Since points in IZ4 always map to IZ3 on the 
next time step, the product M_rMu will always come in 
pairs in the full product for M T . Because this product 
M_rMu produces a diagonal matrix, —I, the full prod- 
uct Eq. ([2]) contains only diagonal matrices, so it is com- 
mutative. This implies that up to a minus sign M r is 
determined by the number of times an orbit visits the 
left half of the unit square, and the parity is determined 
by whether j is even or odd. Putting all of this together 
gives the period r Jacobian stability matrix of a partic- 
ular periodic orbit as 

_ /2™(-l)^ \ 

1V1 r — I q 2~ m (-l) J + 1 / ^ ' 
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with eigenvalues 

|Ac|=2" 
|A e | = 2 ?i 



contracting 
expanding. 



(14) 



The explicit form of the inverse determinant that arises 
in semiclassical calculations is 
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|Det(M T 



(2 m + 2" m -2)- 1 , j odd 
(2 m + 2~ m + 2)- 1 , j even 



(15) 



which may be written as with a ± for notational conve- 
nience, knowing that the sign in front of the 2 is deter- 
mined by the parity of j. Since the main interest is in 
asymptotic calculations (large r) it is sufficient to keep 
the first two terms in the binomial expansion 



(2 L 



± 2)" 1 ps 2 _m - (-l)J'2 -2m+1 



(16) 



It is shown in Sect. IIV Al that actually only the first 
term in the expansion is necessary to investigate certain 
asymptotic fluctuation properties, leaving the approxi- 
mation 
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|Dct(M T - 1)| 



(17) 



Thus, for large period, to leading order |Det(M r — 1)| 
coincides with the contracting stability eigenvalue |A C | of 
the stability matrix. A final note, although the determi- 
nant is the quantity which arises in semiclassical theory, 
some of the classical dynamical systems literature defines 
the uniformity principle with respect to the inverse of the 
stretching exponential fl(| . in which case the correction 
term of Eq. (|16|l is not relevant. 



B. Enumerating the periodic points 

It is quite valuable to be able to count the periodic 
points of fixed period in a given subregion of the unit 
square. To do so we proceed as follows: divide the unit 
square into a grid of boxes of area 2~ fe x 2~ fe , whose 
lower left corners are specified by q = .x\X2--Xk and 
p = .yvyi-.-yk- This is a binary expansion, so each Xi and 
Ui is cither or 1. This is a total of 4 fc boxes. Mem- 
bership of an orbit in a box with a specified lower left 
corner is simply that the orbit has the same first k bits 
for p and q in its binary representation as the lower left 
corner, and arbitrary bits beyond the kth. To start, con- 
sider boxes in the lower right subsquare IZ^. The same 
counting results will hold for regions IZ3 and IZ2 since 
these are merely rotations of 7Z4. The slightly more de- 
tailed counting arguments for region 1Z\ can be found in 
Appendix 

As in Sect. IIII Al the periodic points are written either 
in the form 7m. 7™ or 7m 7m - 7m 7m 1 tne former if j, the 
number of 0-1 or 1-0 transitions in j m , is odd, and the 
latter if j is even. The requirement on the string 7 m for 
a point to be in 71^ is that the first bit is 1. If the last 



bit is 0, then the first form represents a point in IZ4, and 
if the last bit is 1, then it is the second form that is a 
point in IZ4. In either case, m is related to the period r 
by Eq. (fTTjl. Note that j < m - 1, so j < t - 2j - 3 and 
3j < t — 3, so that j is at most L( T /3)J — 1 where |_" • J 
denotes the floor function. The first step is to count how 
many T-periodic points with a fixed value of j there are 
in the box whose corner is specified by q — .\x2--.Xk and 
p = .0j/2— Uk (note the 1 and the are forced because 
the point must lie in IZ4). This point is also represented 
by combining the binary expansions into one expression 
ykyk-i---U2^- 1x2-- -Xk, and similarly for other points. 

For a r-periodic point of the first form, 7 m must look 
like \x2---XkW\...Wiyk---y2^ for some i,k with the restric- 
tion i = T—2j—2k—2 and the total number of transitions 
in this string is j. If the periodic point is of the second 
form, then 7 m = Ix2---Xkw 1 ...w i yk...y2^- 

Let s equal the number of 0-1 or 1-0 transitions in 
Ix2---Xk plus the number of 0-1 or 1-0 transitions in 
0?/2---2/fc- That is, s is the total number of changes for 
the lower left corner point of the box. If s and j are 
both odd, a periodic point would be of the first form, 
and 7m = ^-X2---XkW\...Wiyk---y2^- There are s transitions 
from 1 to Xk and from yk to combined, so there must be 
j — s transitions in the z+1 possible places in w±. ..i/;,?//.. 
i + _ ( t - 2j - 2k - 1 
3 ~s J ~ \ j - s 
this. In the other cases in which j and s may be either 
even or odd, the same result holds and thus we have that 
for each possible value of j, the number of r-pcriodic 
points in a 2~ k x 2~ fe box in H4 with corner value spec- 
r - 2j - 2k - 1 
j - s 

value of j is s, which occurs when all the w's are the 
same as Xk, and the maximum attainable value of j is 
[(r — 2k + s — 1) /3J . In addition, for a given value of s, 
/ 2k — 2 \ 

there are I j possible 2~ k x 2~ fe boxes with s the 

number of transitions in the first k p-bits plus transitions 
in the first k q-bits of the corner point. 

To summarize, the counting just given is for r-periodic 
points in a binary grid of boxes within 7^.4, 7^.3, and IZ2 
where k is the number of bits specifying a box side, and 
s is the number of transitions in the k bits of the q- 
coordinatc of the corner of the box plus the number of 
transitions in the k bits of the p-coordinate of the corner 
point. The index j ranges from s to \_(t — 2k + s — l)/3j 
and for a given j the number of periodic points with that 

r-2j -2k-l 
j-s 

number of period-r points in this box is 



So there are 



ways to do 



ified by s is 



The smallest possible 



value of j is given by 



so the total 



N T ( Sl k) = 



L(T-2fc+s-l)/3j 

E 

3=s 



2j 
j 



2k 

s 



1 



(18) 



Note that for fixed period and box size, the statistics of 
the periodic points within a box are determined entirely 
by the value s of its corner point. Any two boxes with the 
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same value of s will have exactly the same distribution, 
and for each s there arc ( ^ ^ | such boxes. Thus, 

V s / 

the Hannay-Ozorio sum, Eq. ([1}, over all periodic points 
within a binary box in IZ^, TZ^, or IZ2 is 



L(r-2fc+s-l)/3J 

F T (s,k)= 



r-2j- 2k- 1 

3 - s 

1 

|Det(M T - 1)| 



(19) 



and the global form of Eq. (JTJ (excluding Hi) by sum- 
ming over all boxes in H4, 1Z 3 , and IZ2 gives 



2k-2 



The relation of the inverse determinant to the period and 
the symbolic representation of an orbit is, from Eq. (|17[) 
and Eq. (fTTj) . given by 



1 



-T+2j+2 



|Det(M T 



(21) 



which provides an explicit summablc expression for look- 
ing at fluctuations in the uniformity principle. 

The counting arguments for the subsquare 7?.i are 
slightly different, but similar in character, to those pre- 
sented here and the details are given in Appendix [X] In 
fact, the resulting equations arc quite close to the ones 
given in this section. 

Figure [2] shows a plot of all periodic points in the unit 
square at r = 24 and 28. This visualization of the struc- 
ture of the periodic points is interesting in its own right, 
as the points appear to have a fractal-like structure to 
them. In fact the checkered pattern created mimics the 
stable and unstable manifolds of the map. 

The symbolic dynamics and the Markov matrices can 
also be used to find the number of periodic orbits as 
well as to study the uniformity principle. Since the SR 
map is simple enough to permit both a combinatorial 
approach as well as a symbolic dynamics one, it is use- 
ful to present both. Given a periodic point of the first 
type in 7^.4 whose j m string has a binary representation 
IX2X3 • • • W1W2 ■ ■ ■ WiUk ■ ■ ■ 2/20, this translates into an or- 
bit which is a repetition of symbol strings of length r = 
2j+2fc+i+2. Of these 3(s+l) + [2(fe-l)-s] = 2s+2fc+l 
are utilized to specify the fixed corner point. Thus there 
are n = r — (2s + 2k + 1) number of possible "free" sym- 
bols, say Si, . . . S n . A little thought shows that the free 
symbol string has to end with A, that is S n = A, and has 
to be prefixed by an A. Thus it can be either of the form 
CBA ■ ■ ■ A or A ■ ■ ■ A. In both these cases the number of 
periodic points is then given by 



N T ( S ,k)=Y,(To)A Sl (To] 



S1S2 



(To)s n . 



(22) 
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FIG. 2: Fixed points of the time r iterated map. The upper 
square is the plot for r = 24 and the lower square is the 
Rl region for r = 28. The density of fixed points at r = 
28 is roughly 4.6 times that of r = 24. Below, expanding 
Rl renders the lower plot's fixed point density similar to the 
upper plot. The similarity of the expanded Rl fixed point's 
structure to the full phase space at an earlier time illustrates 
the fractal-like structure mentioned in the text. 



N T (s,k) = (T^~ 2k - 2s - 1 ] 



AA ' 



(23) 



Recall that To is the transition matrix in Eq. ([8]). Thus 
the combinatorial problem can be reduced to that of find- 
ing powers of a matrix. From this point the mathematical 
complexity is comparable as both lead to the analysis of 
cubic equations; see Appendix [B] for the combinatorial 
case. 

The uniformity principle sum for the local area can 
also be written compactly in terms of a matrix power, 
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this time the Markov matrix T\. A similar reasoning as 
above leads to 



1 



-2fc-2s-l> 



' AA 



(24) 



Here however, the approximation in Eq. (|21[) is already 
used, as otherwise such a compact formula is not possible. 
Note that with this approximation the global (over the 
whole phase space) uniformity principle sum is simply 
the trace of the power of the Markov matrix. That is 



(25) 



This follows from the fact that the entries in the stochas- 
tic matrix T\ are precisely the multipliers which are ei- 
ther 1/2 or 1. That the stochastic matrix has necessarily 
an eigenvalue 1, and therefore F T — > 1 as r — >• oo is an 
alternative formulation of the uniformity principle. The 
other eigenvalues of T\ whose eigenvalues are less than 
1 in modulus determine both the rate of decay of corre- 
lations as well as approach to uniformity of the periodic 
orbits. We will expand on this below shortly. Almost 
all of the analysis below follows the consequences of the 
binary representation and the combinatorial approach as 
detailed statistics is more transparently done this way. 



IV. STATISTICAL RESULTS 

The results of the previous section and Appendix [B] 
can be used to evaluate sum rule fluctuations. Con- 
sider the local density of the inverse determinant 
|Det(M r - 1)| , which occurs as the natural weight- 
ing for periodic orbits in many semiclassical expressions. 
The first goal is to derive an asymptotic formula for 
its variance. This analysis leads naturally to discussing 
the density and convergence of the remaining component 
F T (s,k) of the Hannay-Ozorio sum rule introduced in 
Sect. Ill B 21 We give an analytic expression for F T (s, k) 
and compute its variance, as well as local and global 
boundaries of convergence for the sum rule. 



A. Local distribution of the inverse determinant 

Consider the regions 1Z±, 7^3, and 7^2 (see Fig.[T]) of the 
unit square whose density of periodic orbits is described 
in Sect. IIII Bl As before, the (similar) discussion for the 
region TZ.\ is left to Appendix [AJ There is a range of 
values taken by |Det(M T — 1)|~ within a local patch of 
phase space, as described in Sect. IIII 51 It was shown that 
the number of period r fixed points with a fixed value of 
the determinant specified by j in a box with parameters 

s and k is given by ^ T ^ ^ ^ where j > s. For 

large period the combinatorial as a function of allowed 
values of j (which can be thought of as a probability 
density for various stability determinant values of fixed 



points in the box) is approximately normally distributed 
(see Appendix [Bj) with an exact mean given by Eq. (|B10I) . 
but written approximately here as 



t -2s -2k -I 
5.148 



s -0.162. 



(26) 



Recall from Eq. (|2lj) that at period r the inverse deter- 
minant may take on the values 2~ T+2 i +2 as j ranges from 
s to [(t — 2k + s — 1)/3J . So in fact, for the discussion of 
|Dct(M r — 1)| , the sum rule F T (s,k) is over terms of 

(t — 2j — 2k — l\ 
the form ( I A J . It turns out that the den- 

V 3- s J 
sity for these terms is normal as before for the t = case; 

note, oddly enough, that does not imply that the density 

for finding a particular value of the inverse determinant 

is lognormal as the convergence with r — 2j — 2k — 1 — > oo 

to normal is too slow. 

Consider the density g(i,t) = ^ n ^ e tl as a func- 
tion of i for a given £, where as in Appendix [B] e* = a. 
Using Stirling's formula to approximate g(i,t), and cal- 
culus, one finds that the maximum value of g(i, t) occurs 
at the value io of i given by 



l 3/3i(i)-2 



(27) 



where /?i(i) is the real root of the cubic equation /3 3 



0. Interestingly this same cubic equation arises 



here for a different problem from the one considered in 
Appendix [B] This method does not give the exact tran- 
sient terms as the recurrence method of Appendix iBldocs. 
but the same structure exists and near the maximum 
at io, g{i,t) is approximated continuously as a Gaussian 
with width on the order of ^fn and so the values of i 
which contribute to the sum are sharply peaked around 
the maximum. Specifically, 



2i 



n+l 



30i 



V27 



(28) 



where of - » (3 ^_ 2)a - 

For the quantity |Det(M T — 1)| _ , for which e* = 4, 
Pi = 2 and so the mean of g(i, t) occurs at io = n/A, even 
though the mean n{n) (for the unweighted combinatorial) 
occurs at about ra/5.148. As n — > oo, the two densities 
tend toward a vanishing overlap since the difference in 
the means grows faster than the widths. In considering 
values taken by the inverse determinant, only those pe- 
riodic points with transition number j which occur near 

t -2s -2k -I 
Jo = - A + s (29) 

contribute to the sum F T (s,k), in spite of the fact that 
there is a vanishing relative fraction of fixed points asso- 
ciated with this value of j, as most points have a value 
of j near (r - 2s - 2k - 1)/5.148 + s - 0.162. 
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B. Important Moments 

We give here explicit expressions for some of the quan- 
tities of interest related to the inverse determinant and 
the SR map, using the results from the methods of Ap- 
pendix [Bj First consider the number of period r fixed 
points within a binary box specified by s and k. The 
number N T (s,k) is given by Eq. (|18j) . which is a special 
case of the sum formula S(n, a) from Appendix [B] with 
a = 1 and n = r — 2k — 2s — l. The form of the solution is 
therefore specified by Eq. (|B7j) . with appropriate values 
for the constants. The topological entropy is given by 
]%t = hi/3i(0), and thus 

N T (s, k) = cie' iT ™ + 2e-' lTn/2 [acos(nd) - 6sin(n0)] 

(30) 

This equation contains the finite-time correction terms to 
the count of fixed points of a binary box, which cannot 
be given by specifying the entropy alone. If r is large, 
the leading term cie hTn dominates and may be used for 
asymptotic calculations. 

The moments for the inverse determinant (which are 
different from the moments involved in the sum rule fluc- 
tuations, ahead) may be computed by averaging over 
powers of |Det(M T — 1)|~ . The most important case 
is the mean, which corresponds to F T (s, k), and gives an 
explicit finite-time correction term to the (infinite time) 
prediction of the uniformity principle. It is also an ingre- 
dient of other sum rule moments. 

From Eq. ([!§)) and Eq. (f2"Tj) . the local form of the sum 
of the inverse determinant over all fixed points within a 
box reads 



F T (s,k) 



L(r-2fc+s-l)/3j 

E 



T~2j - 2k~ 1 \ 2 - T +2]+2 

3 - s ) 

(31) 

which is predicted by the Hannay-Ozorio sum rule to 
asymptotically approach the area of the box, 4~ fe . A 
slight change of variables puts the sum in the generic 
form, Eq. (fBTj) , with a = 4 and n = r - 2k - 2s - 1. The 
solution is thus again of the form of Eq. (|B7[) , where the 
real root /3i of the cubic is exactly 2. After some algebraic 
manipulation, it is seen that F T (s, k) may be written in a 
form which displays both its exponential dependence on 
the period as well as its relation to the phase space area 
as 



F T (s,k) =4" 



acos(n9) — &sin(r 



(32) 

where is the real part of the leading Pollicot-Ruellc 
resonance (here also equal to the positive Lyapunov ex- 
ponent) whose value is given in Eq. (|10[) . The numerical 
values of the constants a, b, and 9 may be calulated using 
the formulas of Appendix. [B] Subtracting the Hannay- 
Ozorio term leaves the oscillating part with time of the 
sum rule as 

F T (s,k) = e- XLT 2 s - k+b/2 [acos(n6) - bsm{n6)]. (33) 



By factoring out the time dependence, it turns out that 
the rate of convergence towards uniformity with increas- 
ing period is exponential, as expected. It is clear from 
the discussion based on symbolic dynamics and Eq. (|2"4")l 
that the rate is governed by the eigenvalues of a Markov 
matrix if it exists or generally by the Pollicott-Ruclle 
resonances. In the case of the SR map, the modulus of 
the eigenvalues of the Markov matrix, which are the res- 
onances, have modulus equal to 1/V2, which is e~ XL . 
It is well-known that the Pollicott-Ruelle resonances are 
generically not related to the Lyapunov exponents and 
therefore the equality for the SR map must be consid- 
ered a coincidence. The close connections between mix- 
ing and uniformity principle makes the emergence of the 
Pollicott-Ruelle resonances as governing the rate of con- 
vergence to uniformity much more natural. 

It is also interesting to consider the convergence bound- 
ary as mentioned in Sect. Ill B 11 This amounts to deter- 
mining the box size (phase space volume) for a given pe- 
riod and location in phase space at which the size of the 
correction term F T (s, k) is just the same order of magni- 
tude as the local area itself. In particular, from Eq. (|33|) . 
if 2 -r/2+s-fc+5/2 = 4 -fe then k = T / 2 - s - 5/2 and the 

volume at the convergence boundary is 



V(s;r)=2 



-r+2s+5 



(34) 



In this way, the local sum rule fluctuations are equally as 
important as the mean, and hence to any results which 
invoke a sum rule on that local scale at that time. Given 
that s varies in the domain 0<s<2fc — 2 or in terms 
of r, < 3s < 2r — 7, the convergence boundary varies 
greatly from one location to another in the phase space, 
i.e. 2~ T < V(s;t) < 2~ T / 3 . Although the local con- 
vergence boundary vanishes everywhere as r — > oo, its 
relative variation tends to infinity. The relatively larger 
boundaries are precisely linked to locally greater inverse 
determinant variation just ahead. 

In Eq. (|16p . the leading correction term to the inverse 
determinant was given, but up to this point not included 
in the calculations. It is important to know if this error is 
subdominant relative to the fluctuating component just 
calculated. If the second expansion term is kept, this 
leads to a sum denoted B T (s, k) of the form 



L(Y-2fc+s-l)/3j 

B T (s,k) = E 



t — 2j — 2k — 1 
x 2- 2r+5 (-16) 3 (35) 



which once again is in the form of the sum discussed in 
Appendix |B] with a = — 16 and n = T — 2s — 2k — 1. 

In fact, more properly, the two most dominant correc- 
tions to the local sum rule are 

F T (s,k) = A~ k +F T {s,k) - B T {s,k) (36) 

We know that the first correction term is governed by 
the Pollicott-Ruelle resonances, but the second term is 
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something else. A priori, it is not obvious which of 
these two correction terms dominates for large period. 
Extracting only the exponential dependence on r gives 
F T (s, k) oc e~ XLT ps (0.71) T . For B T (s, k), because a < 0, 
the dominant fluctuation term comes from the oscillatory 
(a//3)"/ 2 , which is approximately B T (s,k) ps (0.67) t . In 
this case, the first correction term eventually dominates 
over the second, and using the approximation of Eq. (|21[) 
is justified. Had the situation turned out the opposite 
way, then the correction term would not have been a 
Pollicott-Ruelle resonance; we are not aware of an argu- 
ment suggesting that this could not have happened and 
thus both sources of corrections must be considered in 
other cases. 

Before continuing with the spatial fluctuations in the 
sum rule itself, consider the variation of the individual 
inverse determinants contributing to each sum. They 
vary wildly from one fixed point to the next and there is 
greater variation in some regions as opposed to others. 
This gives an s-dependence to their variation within any 
single box. This can be seen by computing the variance. 
The sum of squares of the inverse determinant, Q T (s,k), 
is given by 



(e AT ,s,fc) 



16C1Q / PlQ 



ClN 



4e h 



Pic, 



-2fc-l 



(41) 

which shows asymptotically how the variance varies with 
s, a local characteristic of a particular region of phase 
space (box). Here /3 1Q /e hT is about 1.98 and /3 K3 /4e /tT 
is about 0.495. When the analogous details are worked 
out for boxes in the region IZi, the variance differs only 

by a constant factor of 16 P~Kf \ ■ 

Note that boxes whose lower-left corner has a small 
number of transitions (small s) have smaller variances, 
as well as more r-pcriodic points, than boxes with large 
s. Precisely as found for the local convergence bound- 
aries, the variation of s for moderately large k leads to 
an enormous difference in the variations within different 
boxes of the same size. Although, the variances vanish in 
the limit of r — > oo, the ratios of the variances from one 
box to another another increase indefinitely as the box 
size shrinks. 



Q T (s,k) = ^2 



1 



^ |Det(M T -l)| 

L(r-2fc+s-l)/3j 

E 



T - 2j - 2k - 1 \ 2 _ 2t +4j + 4 

j-s J 

(37) 



This is a sum of the form of Eq. (|B1[) with a = 16 and 
n = t — 2k — 2s — 1. Keeping only the leading term in 
Eq. (|B7|) . the solution is 

Q T {8, k) = A~ T+2s+2 c lQ Pl Q 2k - 2s - 1 (38) 

where /?iq, the real root of the cubic polynomial for 
a = 16, is approximately 2.901 and ciq ps 0.433. The 
subscripts on the numerical constant c and (3 are used 
here to make it clear that these refer to the case Q, for 
which a = 16. The same result can also be derived from 
the symbolic dynamics as 



Q T {s,k) = 



1 



42fe-l \ T 2 



T-2k-2s-l\ 



1 AA 



(39) 



In the notation of Eq. ([5]), the variance of the in- 
verse determinants within a given box is cr 2 (e AT ) = 
M- 2 (e XT ,s,k) - M-i{e XT ,s,k) 2 and 



2 I Ai 

a e 



QrN T 



)-4fc 



N2 



(40) 



C. Sum rule fluctuations 

Next the global variance of the local sum rule F T (s, k) 
is considered due to spatial variation. First we comment 
on the form of the density of F T (s, k). Recall that with 
the method of subdividing the phase space into a grid of 
binary boxes, the value of F T (s,k) locally within a box 
is specified by a parameter s which counts the number of 
to 1 changes in the binary representation of the lower 
left corner of the box (see Sect. IIIIB[) . Thus, F T (s 7 k) 
depends exponentially on s, Eq. (|33|). Furthermore, the 
number of boxes throughout a quarter region of the unit 

f 2k — 2^ 

phase space square with a given value of s is 



as s ranges from to 2k — 2. Thus, the density of the log- 
arithm of F T (s, k) follows a binomial centered at k— 1. As 
with the inverse determinant, the distribution of F T (s, k) 
is described by the product of an exponential function (of 
s, here) and a combinatorial coefficient that is approxi- 
mately normal. A qualitatively similar behavior to the 
discussion of Sect. IIV AI arises in describing the density 
of values taken by the sum formula F T (s, k). 

The variance of F T (s, k) is the average square deviation 
from the mean summing over all boxes and dividing by 
their total number, 4 fc . This is essentially the second 
moment defined in Sect. Ill B 21 



s=0 



2fc-2 
s 



F T ( S ,k) 2 



(42) 



where the leading order of the mean is sufficient. The 
product Q T N T depends on r by a factor (e hT f3\Q / 4) T , 
which is greater than unity. Thus, this term diverges as 
r — > oo. Asymptotically the local variance is just Q T /N T , 



For calculational convenience, Eq. (|33[) is rewritten in the 
form 



F T (s,k) = AeJ s + (Ae^y 



(43) 
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where 7 = In 2 - 291 and A = C2 2~ T / 2 ~ k+3 / 2 e ie{ - T - 2k ~ 1 \ 
Recall that the constants c 2 and 9 arise from the solu- 
tions of the sum formula in Appendix [B] in this case for 
a = 4. This result holds for the regions 724, 723, and 
72-2 of the unit square. For 72i, the expression used for 
F T (s, k) differ only by a constant factor, as shown in Ap- 
pendix [XJ and this factor is accounted for below in giving 
the variance over the entire unit square. 

For large k, it is possible to find a simplified asymptotic 
expression for the variance. Let c 2 = |c 2 |e 4 '' and rj = 
9{t - 2k - 1) + ( giving 



The variance thus decreases exponentially with time, 
again governed by the Pollicott-Ruelle resonance. It also 
increases with the decreasing local volumes. This gives a 
global convergence boundary for the sum rule on which 
the variance over the entire phase space remains a con- 
stant (rather than vanishing). From Eq. (|49| this would 
be given approximately by k 



[Ae^' + iAe^)*} 2 = |c 2 | 2 2^ 2fc+3 [2 2s+1 - 
The expression for the variance becomes 

T2fe-2 



M 2 (F T ,k) = 4- k \c 2 \ 2 2- 



-T-2fc+4 



2fc-2 



-Re(e 21 " 



E 

2k- 



-2Rc(e 



2fc-2 



2ir) 27s 



)]■ 

(44) 



(45) 



Recalling the binomial theorem, each term above may be 
summed explicitly to give 



M 2 (F T ,k) =\c 2 \ z 2- T \Q 



2o-Tic-fc+l [52fc-2 



Letting 1 + e 27 
8cos(4#) gives 



+Rc (e 2 "' [l + e 
1 + 4e- l4e = pe L 



2t 



2fc-2 



(46) 



where p 2 = 17 



M 2 (F T ,k) = \c 2 \ 2 2~ T ~ Ak+A [5 2k - 2 



p 2 cos(2r] + 2u[k - 1])] 



fe-i 



fc-i 



cos(2?7 + 2oj[k - 1]) 



(47) 



Since 49 is not a multiple of 2ir, cos(4(9) is less than unity 
and so is p 2 /25. Thus, the oscillatory terms are subdom- 
inant as k increases. For large k, or small local volume, 
the expression for the variance in each of regions 72-4 , 72-3 , 
and 72-2 becomes: 



M 2 {F T ,k) ~ (3/4)|c 2 | 2 e- 2AiT (5/4) 2fc 



\2fc-2 



(48) 



For the region 72 1 the same equation for F T (s,k) as 
Eq. (|4"5)) applies except that the coefficient A is replaced 
by A' = C2 2- T / 2 -k+i/2 e ze(T-2k+i)_ FrQm tMs it follows 

that the contribution to the variance from the region 1Z\ 
is simply one fourth the value for 72.4, 72-3, or 72-2. The 
asymptotic formula for the variance of F T taken over the 
entire unit square is 



M 2 (F T ,k)~ H| C2 



2 e -2X L r 



2k-2 



(49) 



In 5 /4 



and 



V(r) 



-Tin 2/ In 5/4 



-3. It 



(50) 



V. CONCLUDING REMARKS 



The convergences and fluctuations of classical sum 
rules are interesting in a multitude of ways. Although, 
their corrections may be exponentially suppressed with 
increasing time, the individual contributions can have a 
diverging variance themselves. Another interesting fea- 
ture of local sum rules, as shown herein, is that certain 
fluctuations can be surprisingly large as the location of 
phase space is varied. Correction terms may be related to 
known properties of the system in more general dynami- 
cal systems, such as the topological entropy, the Pollicott- 
Ruelle resonances, or the Lyapunov exponent depending 
on the precise sum rule of interest. It would appear that 
a fluctuation quantity, which depends sensitively on some 
higher power of the stability determinant, if such a quan- 
tity exists, may be more likely to reflect the kinds of 
fluctuations that have been described here on a theoret- 
ical basis for the SR map. The results, however, may be 
suggestive of the type of behavior one might expect in 
a regime where sum rule fluctuations could arise. The 
asymptotic form of several different fluctuation measures 
derived for the Hannay-Ozorio sum, as well as their time 
and length scales, came from the solution of the same 
simple cubic polynomial. The origin of this lies in the 
symbolic dynamics which is a subshift of finite type on 
three symbols and thus there is an equivalent three state 
Markov chain. It may be the case that similar methods 
could be applied to other relatively simple systems, and 
that sum rule corrections could also be derived for these 
systems from the basis of their dynamics. 

The main results for the SR map begin with the calcu- 
lation of the two sources of fluctuations in the Hannay- 
Ozorio sum rule. The first source is governed by the 
dominant Pollicott-Ruelle resonances. It arises from the 
non-uniformity of the locations of fixed points and their 
non-uniform weighting by the leading behavior of their 
inverse stability determinants. The second source arises 
from the effects of next-to-leading order corrections to 
the inverse stability determinants. These corrections are 
not governed by the Pollicott-Ruelle resonances, but are 
also exponentially decreasing in time. The dominant cor- 
rection here comes from the first source and hence the 
Pollicott-Ruelle resonances, however we do not currently 
know whether this must be the case for general chaotic 
dynamical systems. 

It is a matter of how closely the stretching multipliers 
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approximate the determinant in Eq. (|15|) . and it could 
be that for some other chaotic system they are different 
enough to produce corrections that dominate the one due 
to the resonances, although these would still be present. 
In the specific case of the SR map the second term in 
Eq. (|16[) . the principal correction, is an oscillating sum 
because the periodic orbits are reflecting hyperbolic if the 
number of — 1,1 — bonds are odd. This term can be 
written as a trace of the power of the matrix 

1/4 1/4 \ 

-10 (51) 
10/ 

which is different from the T 2 matrix in that the element 
(2, 1) is — 1 rather than 1. This ensures that each time 
the orbit gets rotated, it acquires a negative sign. Al- 
ternately, for each CBA part of the symbolic string a 
negative sign is acquired. 

The leading eigenvalue of this matrix has a modulus of 
ss 0.67 which is smaller than the subleading eigenvalue, 
« 0.707, of 2~i that gives the Pollicot-Ruelle resonance. 
If the orbit were all direct hyperbolic then in the above 
matrix the —1 will change to 1 and this will be same 
as T2 whose leading eigenvalue is 0.725 which is larger 
than 0.707, and would have dominated the corrections. 
The fact that some of the orbits are reflecting hyperbolic 
seems to have been crucial to lower the contribution from 
corrections that come from the fact that a det (J — I) is 
present instead of just the multipliers. 

However if the sum rule is weighted by (the inverse) 
of the largest eigenvalue of the stability matrix, the 
Pollicott-Ruelle resonances will govern the corrections 
to the sum rule, especially if there is a finite symbolic 
dynamics description of the system. The relevance of 
Markov and related matrices (To, T\, T2, . . .) for the cal- 
culation of the fluctuations indicate possible connections 
with the Thermodynamic Formalism especially as ap- 
plied to finite Markov processes [l7| ■ 

A second result shows how the relative local variations 
of the inverse determinants varies infinitely broadly at 
long times. Finally, the relative variation of the sum rule 
applied locally also has an infinite width while maintain- 
ing an exponential convergence rate for fixed phase space 
volume. We gave convergence boundaries that show how 
small a local volume may be considered for a given time of 
propagation if one expects convergence to the asymptotic 
sum rule result. Again, the relative size of a converged 
local volume depended on location and varied infinitely 
broadly while maintaining exponential convergence with 
time at fixed volume 

It would be extremely interesting to investigate other 
sum rules, especially those that connect to quantum fluc- 
tuation properties of eigenfunctions and transport. The 
various localizing effects giving rise to eigenfunction scar- 
ring [l8j , localization manifestations of time scales intro- 
duced by transport barriers |19|. and interaction effects 
linked to Friedel oscillations j20L [2lj give a few inter- 
esting directions for further studies. As mentioned ear- 



lier, the SR map is easily studied quantum mechanically 
and would be one possible way to study sum rules aris- 
ing from quantum fluctuations properties involving eigen- 
functions. 
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Appendix A: Region IZi 

Here the counting arguments and several results for the 
region IZi of the unit square, which has mostly been ig- 
nored in the body of the text, are presented. The reason 
for leaving this discussion here is that many of the de- 
rived results closely resemble those for the other regions, 
although the arguments are somewhat longer. 

We begin with an extension of Sect. IIII Bl bv counting 
the number of period r points in a binary 2~ k by 2~ k 
box in the region 7£i (Fig. [IJ, but not on the bottom 
row of boxes. The lower left corner of each box is 
defined by 0yk---y20-0x2---XkQ, with the condition that 

not all of the y's are zero. The upper right corner is 
given by lj/fc...y20.0x2...Xfe 1. It is easy to see that after 

applying the inverse transformation of the map some 
number of times, each square of area 4~ fc will be mapped 
into a rectangle in 7Z4 of the same area, although not 
square, and also that the upper right corner of the 
square m 1Z\ gets mapped into the lower left corner 
of the rectangle in IZ4. Thus the box in IZi with 
lower-left corner ?/&... 10. ..0.0cc2....Tfc (where not all the y's 
are zero) has upper right corner lyfc...lO...0 . 0x2---Xk 1 

which transforms under the inverse transformation as 
follows: lj/fc...lO...0 . 0x2---Xkl -s— lj/fe...l.0... 0x2---Xkl <— 
Qx k ...x 2 l...l .l...y k l <~ Q£/fe...O . l...lx2...SfcO, so 
the lower left corner of the rectangle in IZ4 is 
2/&...0 . 1...lx2---Xki with transition number one less 
than the lower left corner of the box in 1Z\ . 

Let s be the number of transitions of the lower left 
corner of a box in 1Z\. Then by the same counting argu- 
ment that was used before (the fact that it is a rectangle 
instead of a square does not change things), the number 
of r-periodic points in this box in 1Z\ with a given tran- 

( r — 2j — 2k — l\ 
sition number 7 is I . , , as 7 ranges from 

V 3-s+l J 

s - 1 to [(r - 2k + s - 2)/3j . So this looks just like it did 
for the IZ4 case except s is replaced by s - 1. 

Now, the lower left corner of a box that we are consid- 
ering in this region has, say, t transitions in the position 
coordinate and v transitions in the momentum coordi- 
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nate, where < t < k — 1, and 1 < v < k — 1, and 
s = t + v. The combinatorial expression above shows that 
the number of r-periodic points in a box depends only on 
the total s and not how it is distributed between t and 
v, however it is necessary to consider t and v separately 
because it is v that is restricted to be greater than zero 
and not just the sum of t and v. We separately choose t 
transitions from k — 1 places for position, and v transi- 
tions from k — 1 places for momentum, with v restricted 
to be greater than zero. So, for example, the local form 
of the sum of the inverse determinant for boxes in IZi 
excluding the bottom row (denoted lZi a ), analogous to 
expression Eq. (fT9|) , is 



L(T-2fe+t+ti-2)/3j 

F T {t + v,k;K la ) = 

j=t+v-l 



t - 2j -2k -I 
j -(t + v-1) 

1 

|Det(M T - 1)| 

(Al) 



and the sum over all boxes in lZi a , analogous to expres- 
sion Eq. ([20]), is 



v=l t=0 



(A2) 

Considering the bottom row of boxes in IZi (denoted IZu) 

a similar, but more tedious argument, gives the number 

( t — 2j — 2k — 1 \ 
of period r points again as I j s i J as this 

time j ranges from s to [(t — 2k + s — 2)/3j . So the local 
sum of the inverse determinant for boxes on the bottom 
row of IZi is 



For the sum of squared deviations over boxes in IZi the 
expression to evaluate is 



v=l t=0 
k-1 



(A6) 



where F r (s) = F T (s) - 4" fc : . 

There is an identity due to Vandermonde [13] which 
simplifies the double sum above and gives a result that 
is almost exactly like the sum for region 7?4. Denoting 
the double sum by W gives 



w-EE(*; 1 )(*7 1 J*<*+<'>" 

v=l t=0 v ' v 



k-1 v+k-1 

E E 

V—l S—V 







1 









Fr{sf (A7) 



where s — t + v. Interchanging the order of summation 
and breaking this into two sums generates 



k-1 



^-E^) 2 E 

s=l v=l 

2k-2 k-1 

E^(^) 2 E 



k-i\ f k-i 



k-i\ f k-i 



s—k v— s — fc+1 

Vandermonde's convolution identity is 



s — v 



L(r-2fc+.s— 2)/3J 

F T { s ,k-n lb )= 



t — 2j — 2k — 1 
j-s + 1 

1 



(A3) 



|Det(M r - 1)| 
and the sum over all boxes on the bottom row of 7li is 

k-1 



iF Tb - 



k-1 

•s 



F T { s ,k-n lb ) 



(A4) 



The more complicated sums that arise for region 1Z\ 
may be simplified for certain calculations of interest. In 
particular, consider the calculation of the variance for the 
Hannay-Ozorio sum in Sect. HVCl Equations ([AT]) and 
(|A3|) . like their 72-4 counterpart Eq. (fT9|) . are sums of the 
form of Eq. (jBlj) from Appendix |B] with a = 4, and they 
may both be expressed as 

F T (s;TZ 1 )^A- k + A'e^ s (A5) 
where A' = C2 2-^/ 2 - fe + 1 /2 e ie(r-2fc+i) and 1 = \ n2 - 29i. 



c \ / a+c 
b I \ d-b I ~ I d 



(A8) 



where the sum is over all values of b for which the sum- 
mand is not zero. This gives 



E 

v=l 



k-1 \ ( k-1 



2k -2 \ f k-1 



when 1 < s < k — 1, where the subtracted term corre- 
sponds to v = 0. Also note that 



k-1 

E 



k-i \ f k-i 



2k -2 
s 



v—s — k-\-l 

when k < s < 2k — 2. Combining the two gives 



2k— 2 * . \ k—1 / \ 

^=E( 2fc ; 2 H(*) 2 -Ef"; 1 )^^ 
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and therefore, 



2fc-2 



a 2 (F T ) = 4- fc [£ ( 2k : 2 )F T (sf 



8 = 1 

k-1 



Ei*: 1 )^) 2 



s=l 

fc-1 



E( fc ; x )^) 2 ] 



or 



2fc-2 



s=l 
2fe-2 



a 2 (F r )=4- fe [£ ( 2fc ; 2 )F T ( S ) 2 + F r (0) 2 i 



which is simply 



s=0 



2fc-2 



4- fe E ( 2fc : 2 iw 



(A10) 



(All) 

This sum is of exactly the same form as Eq. ((32) for 
computing the variance for the other regions of the unit 
square. 



Appendix B: A sum formula for the SR map 

Upon examination of the form of Eq. (|T9|) and given 
the result of Eq. ([21]) , it happens that in order to arrive 
at closed form expressions for fluctuations in the Hannay- 
Ozorio sum, Eq. ([T]), it turns out that several sums of the 
form 



L«/3J 

S(n,a) = ^2 



i=0 



2i 



(Bl) 



are needed for various real values of a, with n a positive 
integer. This section gives a general discussion of such 
sums and presents a method of finding closed form ex- 
pressions for them before moving on to the main results 
of interest. In the analysis of the SR map, the cases a = 
1, 4, 16 and -16 show up naturally when considering the 
lower order moments of Sect. Ill B 21 

Obtaining a closed form for the sum may begin by 
finding a recursion formula for it, and then using a 
standard technique for solving such recursions. Recall 

f k 

the recursion for building Pascal's triangle: 



fc-1 
i - 1 



fc-1 



where fc and i are greater than 



one. Applying this gives 

L«/3J 



S(n, «) = E 1 



i=0 
|n/3J 

E 

i=l 



n — 2i 
i 



n-H-\ 

i - 1 



L(n-X)/3J 



+ E 

i=0 
L(n-3)/3j 

- E 



n-2i-\ 



n - 3 - 2j 



3=0 

L(n-X)/3J 



E 



n - 1 - 2i 



a 1 (B2) 



where j = i — 1 in the first summation. This produces 
the recursion relation 

5(n, a) = S(n — 1, a) + aS(n — 3, a) (B3) 

for n > 3 with initial conditions S(0,a) = S(l,a) = 
5(2, a) = 1. Of historical note, the 14*' 1 century In- 
dian mathematician Narayana studied a problem of the 
proliferation of cows (each offspring gives birth after its 
third year) that leads to this very same recursion relation 
with a = 1 [23|. A standard technique for solving such 
recurrence relations is to look for solutions of the form 
5(n, a) = /3 n , and thus 5(n-l, a) = /T 1 /?", 5(n-3, a) = 
f3~ 3 /3 n . Plugging these expressions into Eq. (|B3|) the fac- 
tor /3 n cancels and leaves the cubic equation 



a = 0. 



(B4) 



The three roots of this cubic /3i, (3%, P3 give three so- 
lutions of the recurrence /3™, j3%, and (3^ . The differ- 



ence equation Eq. (|B3[) is third order, linear, and ho- 
mogeneous, and the standard theory of such difference 
equations (analogous to that for differential equations) 
says that if there are three linearly independent solutions, 
then the general solution may be formed as a linear com- 
bination of these independent solutions. Naturally this 
cubic equation also appears when using the matrices from 
symbolic dynamics. Indeed the characteristic equations 
for To, T\ and T2 are, up to a scaling, the same as the 
cubic equations with a = 1, 4 and 16, respectively. 

The cubic polynomial /(/?) = (3 3 — 2 — a has a local 
maximum of —a when (3 = and a local minimum at 
/3 = 2/3. It has one real root, say fix, when a > and 
also when a < —4/27, which covers all of the cases of 
interest. The other two roots are complex conjugates, 
/?2 = re l ° = S + and /?3 = re~ %B = S — ij. Since 
(3f((3i — 1) = a, it implies that f3i > 1 when a > 0, and 
f3i < —1/3 when a < —4/27, so that (3\ has the same 
sign as a. 

Because the three roots are distinct, the three solutions 
are independent and the general solution to Eq. (|B3[) may 
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be written as 



which gives a complete and explicit solution to Eq. (|B1|| 



S(n,a) = dft + c 2 {re lS ) n + c 3 {re~ te ) r 



(B5) 



where the coefficients c\, c 2 , C3 may be found from the 
initial conditions on the recurrence (in all cases here the 
initial conditions on S(n, a) are real). 

It is possible to express the two complex roots, as well 
as the coefficients, in terms of the real root Pi and in 
terms of a. The constant c\ is real and C3 is the conjugate 
of C2, which can be denoted as c 2 = a + ib, C3 = a - ib 
where a and b are real. The polynomial of Eq. (|B4[) may 
be factored as (/? — Pi)(P — S — ij)(P — S + 17). 

Comparing the constant terms of the polynomial writ- 



ten both ways gives Pi(5 + 7 ) 



Thus, the magni- 



tude of the complex roots is r = (£ 2 +7 2 ) 1 / 2 = {a/ Pi) 1 / 2 . 
Because P\ — /3± = a/ Pi, it is also true that r < fix when 
a > 0, but r > f3i when a < —4/27. The significance is 
that for large n the oscillatory terms in Eq. (|B5[) are dom- 
inated by the first term when a > 0, but the oscillatory 
terms are dominant when a < —4/27. This observation 
is used ahead in deriving several asymptotic results. 

Comparing the coefficients of the square terms gives 
5 = (1 — /3i)/2, which is a negative number when a > 
and a positive number when a < —4/27. The two com- 
plex roots are in the second and third quadrants when 
a > 0, and in the first and fourth quadrants in the other 
case. For a > 0, we can take 8 = it + arcsin((l — (3i)/2r) 
which lies in the second quadrant, and for a < —4/27 
we can take 8 = arcsin((l — /?i)/2r) which is in the first 
quadrant. 

With the complex roots in terms of the one real root, 
it suffices to find the real root, which can be expressed 
straightforwardly for the regime of interest, i.e. either 
a > or a < -4/27. In that case, with x = 1 + 27a/2, 



1/3 



1 + (x+ \Jx 2 - 1 ) + [x - \Jx 2 - 1 



1/3' 



A-5 

(B6) 

Rewriting Eq. (|B5| in terms of the real constants ci, a, 
b gives 



The counting results for periodic points in the previous 
section beg the question, 'what is the asymptotic density 

f n — 2i \ 

of the combinatorial expression I . I as a function 

of iV . It turns out that it is possible to find the asymp- 
totic mean, variance, and density of i (as n approaches 
infinity) by essentially the same algebraic methods used 
in the recurrence relation. Using the moment-generating 
function, or by using Stirling's approximation (essentially 
a saddle point expression), it can be shown that the den- 
sity of ( U ^ ^ ^ , when properly normalized, converges 

to a normal density. The moment-generating function 
technique also gives simple formulas for the asymptotic 
mean and variance. More specifically, 



S(n,e f ) 



Ln/3j 

E 

i=0 



2i 



(B9) 



where one substitutes a = e , and the moment- 
generating function for this density is 4>(n, t) = 
S(n, e*)/5i(n, 0). Large n gives S(n, e*) -> c(t)pi(t) n and 
(f>(n,t) = c(t)(3i(t) n [c(0)/3i(0)™] _1 where is the real 
root of P(t) 3 - (3(t) 2 - e* = 0. The details are omitted, 
but by differentiating the cubic equation all of the deriva- 
tives of 4>(n, t) can be found, which can be used to find 
the moments of the density. When a = e* = 1 the real 
root of Eq. (|B4j) is /3i(0) = 1.46557123 . . ., which in the 
next section is seen to have special significance to this 
map. The asymptotic mean of this density can be shown 
to be 



3 + 01(0) [ 



n-2 



3 + /? 2 (0) 



(B10) 



f a\ n/2 

S(n,a) = cipl + 2 i-^j [acos(n6») - bsm(nd)} 

(B7) 

Putting in the initial conditions gives three real equations 
for the coefficients which may be solved in terms of a and 
Pi. Skipping the algebraic steps, one finds 



ci 



a + Pl 
3a + PI 

Pi 



3a + P( 



(3a + p() \3Pi + 1 



1/2 



(B8) 



and the variance 



a 2 (n) 



[3 + ft 2 (0)] 
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(0) 



(Bll) 



The scaling of the mean is n and the width is n 1 / 2 . All the 
higher reduced cumulants (rescaled by the appropriate 
power of the width) vanish in the limit of n — > 00. 
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